In [7]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

In [8]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

In [9]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

In [10]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

In [11]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

In [12]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

In [13]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

In [14]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

In [15]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

In [16]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

In [17]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab05.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab05.py:44: RuntimeWarning: divide by zero encountered in double_scalars

turning = (v0/v_list[-2] < 0.0)

In [18]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 10/E21.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 10')

2.5025864848083885 2.5036274076761975

/Users/amandaash/Desktop/PHYS_3210/Week 10/E21.py:106: UserWarning: Matplotlib is currently using module://ipykernel.pylab.backend_inline, which is a non-GUI backend, so cannot show the figure.

fig1.show()

/Users/amandaash/Desktop/PHYS_3210/Week 10/E21.py:109: RuntimeWarning: divide by zero encountered in log10

stability = -np.log10(np.abs((E_tot-E_tot[0])/E_tot[0]))

1.6700847418249443 1.674165797300528

/Users/amandaash/Desktop/PHYS_3210/Week 10/E21.py:106: UserWarning: Matplotlib is currently using module://ipykernel.pylab.backend_inline, which is a non-GUI backend, so cannot show the figure.

fig1.show()

/Users/amandaash/Desktop/PHYS_3210/Week 10/E21.py:109: RuntimeWarning: divide by zero encountered in log10

stability = -np.log10(np.abs((E_tot-E_tot[0])/E_tot[0]))

1.2555006764821623 1.2598267915100132

/Users/amandaash/Desktop/PHYS_3210/Week 10/E21.py:106: UserWarning: Matplotlib is currently using module://ipykernel.pylab.backend_inline, which is a non-GUI backend, so cannot show the figure.

fig1.show()

/Users/amandaash/Desktop/PHYS_3210/Week 10/E21.py:109: RuntimeWarning: divide by zero encountered in log10

stability = -np.log10(np.abs((E_tot-E_tot[0])/E_tot[0]))

1.007824640651541 1.009577543176985

/Users/amandaash/Desktop/PHYS_3210/Week 10/E21.py:106: UserWarning: Matplotlib is currently using module://ipykernel.pylab.backend_inline, which is a non-GUI backend, so cannot show the figure.

fig1.show()

/Users/amandaash/Desktop/PHYS_3210/Week 10/E21.py:109: RuntimeWarning: divide by zero encountered in log10

stability = -np.log10(np.abs((E_tot-E_tot[0])/E_tot[0]))

0.8429992525831275 0.8426868158909813

/Users/amandaash/Desktop/PHYS_3210/Week 10/E21.py:106: UserWarning: Matplotlib is currently using module://ipykernel.pylab.backend_inline, which is a non-GUI backend, so cannot show the figure.

fig1.show()

/Users/amandaash/Desktop/PHYS_3210/Week 10/E21.py:109: RuntimeWarning: divide by zero encountered in log10

stability = -np.log10(np.abs((E_tot-E_tot[0])/E_tot[0]))

0.7254922074480112 0.7236414065808147

/Users/amandaash/Desktop/PHYS_3210/Week 10/E21.py:106: UserWarning: Matplotlib is currently using module://ipykernel.pylab.backend_inline, which is a non-GUI backend, so cannot show the figure.

fig1.show()

/Users/amandaash/Desktop/PHYS_3210/Week 10/E21.py:109: RuntimeWarning: divide by zero encountered in log10

stability = -np.log10(np.abs((E_tot-E_tot[0])/E_tot[0]))

In [19]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 0.683730225975385 N

In [20]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 2.6810623217992444 N

In [21]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 2.6810623217992444 N

In [22]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 2.6810623217992444 N

In [23]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 2.6810623217992444 N

In [24]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 2.6810623217992444 N

In [25]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

In [26]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

In [27]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

In [28]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

In [29]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

In [30]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

In [31]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

In [32]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

Traceback (most recent call last):

File "<ipython-input-32-43bbbc208fdb>", line 1, in <module>

runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

File "/Users/amandaash/anaconda3/lib/python3.7/site-packages/spyder_kernels/customize/spydercustomize.py", line 827, in runfile

execfile(filename, namespace)

File "/Users/amandaash/anaconda3/lib/python3.7/site-packages/spyder_kernels/customize/spydercustomize.py", line 110, in execfile

exec(compile(f.read(), filename, 'exec'), namespace)

File "/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py", line 161, in <module>

natural_period, natural_frequency = find_period(v_natural)

TypeError: find_period() missing 1 required positional argument: 'time_step'

In [33]:

In [33]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Reloaded modules: __mp_main__

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

Traceback (most recent call last):

File "<ipython-input-33-43bbbc208fdb>", line 1, in <module>

runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

File "/Users/amandaash/anaconda3/lib/python3.7/site-packages/spyder_kernels/customize/spydercustomize.py", line 827, in runfile

execfile(filename, namespace)

File "/Users/amandaash/anaconda3/lib/python3.7/site-packages/spyder_kernels/customize/spydercustomize.py", line 110, in execfile

exec(compile(f.read(), filename, 'exec'), namespace)

File "/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py", line 161, in <module>

natural_period, natural_frequency = find_period(v_natural)

TypeError: find_period() missing 1 required positional argument: 'time_step'

In [34]:

In [34]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

2.679168218991807

In [35]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

2.679168218991807

In [36]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

2.679168218991807

In [37]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

2.679168218991807

In [38]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

2.679168218991807

In [39]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

2.679168218991807

In [40]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

3.162464922075491

In [41]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

3.162464922075491

In [42]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

2.679168218991807

In [43]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

Traceback (most recent call last):

File "<ipython-input-43-43bbbc208fdb>", line 1, in <module>

runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

File "/Users/amandaash/anaconda3/lib/python3.7/site-packages/spyder_kernels/customize/spydercustomize.py", line 827, in runfile

execfile(filename, namespace)

File "/Users/amandaash/anaconda3/lib/python3.7/site-packages/spyder_kernels/customize/spydercustomize.py", line 110, in execfile

exec(compile(f.read(), filename, 'exec'), namespace)

File "/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py", line 158, in <module>

x_natural, v_natural, t_natural = harmonic_oscillator_friction_beta(p_value,k_constant,v_initial,x_initial,mass,dt,t_initial,t_final,0,0,0)

File "/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py", line 24, in harmonic_oscillator_friction_beta

if np.isclose(v,0, atol = 10**-5):

File "/Users/amandaash/anaconda3/lib/python3.7/site-packages/numpy/core/numeric.py", line 2523, in isclose

if all(xfin) and all(yfin):

File "/Users/amandaash/anaconda3/lib/python3.7/site-packages/numpy/core/fromnumeric.py", line 2245, in all

return _wrapreduction(a, np.logical_and, 'all', axis, None, out, keepdims=keepdims)

File "/Users/amandaash/anaconda3/lib/python3.7/site-packages/numpy/core/fromnumeric.py", line 84, in _wrapreduction

return reduction(axis=axis, out=out, **passkwargs)

File "/Users/amandaash/anaconda3/lib/python3.7/site-packages/numpy/core/_methods.py", line 46, in _all

return umr_all(a, axis, dtype, out, keepdims)

KeyboardInterrupt

In [44]:

In [44]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

3.162464922075491

In [45]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

3.162464922075491

In [46]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

Static Friction failure, F_f = 4.4145 N, F_spring = 1.1127872145355915 N

3.162464922075491

In [47]: runfile('/Users/amandaash/Desktop/PHYS_3210/Week 11/Lab_5.py', wdir='/Users/amandaash/Desktop/PHYS_3210/Week 11')

3.162464922075491

In [48]: 0.1/(2*np.pi)